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The variable solar magnetic activity known as the 11-year solar cycle has the longest history of solar observations. These cycles dramatically affect conditions in the heliosphere and the Earth’s space environment. Our current understanding of the physical processes that make up global 
solar dynamics and the dynamo that generates the magnetic fields is sketchy, resulting in unrealistic descriptions in theoretical and numerical models of the solar cycles. The absence of long-term observations of solar interior dynamics and photospheric magnetic fields hinders 
development of accurate dynamo models and their calibration. In such situations, mathematical data assimilation methods provide an optimal approach for combining the available observational data and their uncertainties with theoretical models in order to estimate the state of the solar 
dynamo and predict future cycles. In this presentation, we will discuss the implementation and performance of an Ensemble Kalman Filter data assimilation method based on the Parker migratory dynamo model, complemented by the equation of magnetic helicity conservation and long- 
term sunspot data series. This approach has allowed us to reproduce the general properties of solar cycles and has already demonstrated a good predictive capability for the current cycle, 24. We will discuss further development of this approach, which includes a more sophisticated 
dynamo model, synoptic magnetogram data, and employs the DART Data Assimilation Research Testbed. 


Components of the Solar Cycle Prediction Approach Preliminary Analysis of Prediction Solar Cycle 25 Uncertainties 


Variations of solar activity are a result of complicated dynamo processes in the convection zone. We consider this | To characterize the level of solar activity during systematic observations of sunspots, the relative sunspot Previous experience and tests have shown that the EnKF procedure 
phenomenon in the context of sunspot number variations, for which we have detailed observational data during the © Number was introduced by Wolf (1850): R = k(10g + n), where k is a correction factor, depending on based on the dynamo model and sunspot number measurements has 
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past 23 solar cycles. However, despite the known general properties of solar cycles, a reliable forecast of the 11- | observing conditions, g is the number of sunspot groups, and n is the number of individual sunspots. good predictive capabilities for estimating future solar activity in a time 
year sunspot number is still a problem. The main reasons for these forecasting uncertainties are imperfect dynamo a range from 7 - 8 years to a whole solar cycle. However, our attempts to 
models and deficiency of the necessary observational data. a) a ee ‘original’ sunspot data series == tC Pep meget ae ea ie make predictions for a period longer than one cycle often fail due to 
Shes Wells Dynamo model jee - v2.0 sunspot data series A accumulation of errors. In this section, | consider in more detail the 
——_ - Parker 1955. Kleeorin & Ruzmaikin. 1982 250- | | < ie possibility of early forecasts of sunspot cycles and present an _ initial 
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Comparison of the annual sunspot number series: the original (marked ‘original’ solid curve), which a ~ — + — 
was used for the solar cycle 24 prediction by (Kitiashvili & Kosovichev 2008), and the ‘revised’ 1980 1990 2000 2010 


: . year 
version (v2.0, Source: WDC-SILSO, Royal Observatory of Belgium, Brussels) by applying a new Synoptic magnetogram. The color scale is saturated at +/-15G. The white 


2 calibration to the historical and modern data sets (dotted curve, Clette et al. (2015) is shown in panel dashed lines indicate different moments of time: 1992: 1997 and 2015. 
5 ie ae a). Panel b) shows the ratio of these data sets. 
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Step 4: Observation increments are computed based on the relationship between 
the observed value and expected error vs. the distribution of expected mei ox | timp_tlast_pbs = 1991 a time_last_ obs = 1992 gol. He test obs 1294 
observation values. 1990 1995 2000 = 2005 | | | 

Step 5: The increments are regressed onto the model states and the state values — 
adjusted. 


— Comparison of the sunspot number prediction for Solar Cycle 24 (red 
Step 6: Each model state is advanced in time and the cycle repeats until curve, Kitiashvili & Kosovichev, 2008) and actual observations of 
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Testing the Prediction Capabilities solution according to annual sunspot number (green diamonds). 
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show the best EnKF estimate of the model variations. Panels d-f and j-/ show _ toroidal 


Prediction of solar cycles is one of most interesting problems closely linked to dynamo processes inside the Sun. Numerous earlier attempts to predict future solar cycles were mostly based magnetic field errors of the model for each data assimilation case. 


on empirical relations derived from observations of previous cycles and provided a wide range of predicted strengths and durations of the cycles. The difficulty is due to our incomplete 
understanding of the physical mechanisms of the solar dynamo and also due to observational limitations that result in significant uncertainties in the initial conditions and model parameters. a Foie abla icantar fis inchicsneh rere eclecicaecen Gaon 
. , . . . : . . naerson J. . Bulletin O e American WieteoroiogiCal SOCciety, ‘ 
We have developed a relatively simple non-linear mean-field dynamo model, which nevertheless can describe essential general properties of the cycles and the observed sunspot number jderson JL, Collins N. 2007. Journal of Atmospheric and Oceanic Technology, 24, 1452 
series (such as Waldmeier’s rule). Combined with the data assimilation approach, this model provides reasonable estimates for the strength of the following solar cycles. In particular, the  Clette, F., Svalgaard, L., Vaquero, J. M., & Cliver, E. W. 2015. Space Sciences Ser. of ISSI, Vol. 53, 35 


prediction of Cycle 24 calculated and published in 2008 is holding quite well so far. Evensen G. 2007. Data Assimilation. Springer | 
eer 238 : . : ; ; Kalman R.E. 1960. Problems. J. Basic. Eng., 82, series D, 35-45 
The initial prediction of Cycle 25 shows that this cycle will start in about 2019 - 2020, reach the maximum in 2023 - 2024, and the mean sunspot number at the maximum will be ~90 (for the — Kitiashvili |. N., Kosovichev A. G. 2008. ApJ, 688, L49-L52 


v2.0 sunspot number series) with an estimated error of ~15%. The test simulation runs for early cycle predictions have been performed and identified the following criteria for a successful Sail a Sacer ae ma eles a aout 93-68 
prediction capability: a) the model errors relative to observations should be less than 20% for the last 10 years, and b) the prediction should be performed starting from a period when either Giscivi IN. 2016. ApJ, 837 7 a eal 
the toroidal or poloidal field is dominant relative to the other. This corresponds to two moments of time: the polar field reversals shortly after the solar maxima (strong toroidal field and weak  Kleeorin N.I., Ruzmaikin A.A. 1982. Magnetohydrodynamics, 18, 116 


poloidal field) and during the solar minima (strongest poloidal and weak toroidal fields). Parker E.N. 1955. ApJ, 122, 293 
oa ve a Schatten K. 2005. Geophysical Research Letters, Vol. 32, CitelD L21106 
The next steps of this study are (1) transition from the sunspot number to physical quantities such as synoptic magnetograms and magnetic helicity and (2) using more advanced dynamo _ Ssvalgaard L., Cliver E.W., Kamide Y. 2005. ASP Conf. Series, Vol. 346, p.401 


models that provide a better physical basis for development of a high-fidelity pipeline for solar activity predictions. Wolf, R. 1850. Astronomische Mitteilungen der Eidgen ossischen Sternwarte Zurich, 1, 15 
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